Aspects of superfluid cold atomic gases in optical lattices 
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We review our studies on Bose and Fermi superfluids of cold atomic gases in optical lattices at 
zero temperature. Especially, we focus on superfluid Fermi gases along the crossover between the 
Bardeen-Cooper-Schrieffer (BCS) and the Bose-Einstein condensate (BEC) states, which enables us 
to study the Bose and Fermi superfluids in a unified point of view. We discuss basic static and long- 
wavelength properties (such as equation of state, incompressibility, and effective mass), energetic 
stability, and energy band structures of the superfluid Fermi gases in an optical lattice periodic along 
one spatial direction. The periodic potential makes pairs of atoms to be strongly bound, and this 
can affect the static and long-wavelength properties and the stability of the superflow. Regarding 
the band structure, a peculiar loop structure called "swallowtail" can appear in superfluid Fermi 
gases as well as in the Bose case while its mechanism of the emergence is very different from the 
bosonic counterpart. Other quantum phases that the cold atomic gases in optical lattices can show 
are also briefly discussed from their role as quantum simulators of the Hubbard models. 
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I. INTRODUCTION 

In 1995, the atomic physics witnessed the coolest state 
of matter realized in the laboratories to that day. With 
lasers and electromagnetic fields, atomic Bose gases were 
trapped and cooled down to within a micro-Kelvin to 
make them coalesce into a single quantum wave, known 
as a Bose-Einstein condensate (BEC) state [1]. With 
the technology and tools developed and upgraded, fi- 
nally about 10 years later, physicists succeeded in cooling 
atomic Fermi gases down to the even lower temperatures 
and in providing a conclusive evidence that the interact- 
ing atomic Fermi gas could be controlled to be in the 
state from Bardeen-Cooper-Schrieffer (BCS) superfluid 
of loosely bound pairs to a BEC of tightly bound dimers 
[2]. This quantum wonderland of cold atomic gases has 
been proven to be rich in physics and it has also shown 
us new realms of research [3-7]. Unprecedented high- 
precision control over the cold atomic gases has made it 
possible to mimic various quantum systems. It has also 
paved a way to new physical parameter regimes which 
had not been attainable in other quantum systems and 
therefore had not been well thought of. 

In the cold Fermi gas with an equal mixture of two hy- 
perfine states (for simplicity, called as spin up/down), 
wide-range control of the short-range interatomic in- 
teraction (characterized by an s-wave scattering length 
fls) via Feshbach resonances [8], has enabled mapping 
a new landscape of superfluidity, known as "BCS-BEC 
crossover", smoothly connecting from the BCS super- 
fluid to the BEC superfluid [9]. For weakly attractive 
interaction (l/fcpas <SC —1, with fcp conventionally being 
the Fermi momentum of the free Fermi gas of the same 
density), the gas shows the superfluid behavior originat- 
ing from the many-body physics of Cooper pairs which 



are formed from weak pairing of atoms of opposite spin 
and momentum and whose spatial extent is larger than 
the interparticle distance. For strongly attractive in- 
teraction {1/k-pas » 1), the gas shows the superfluid 
behavior which can be explained by the Bose-Einstein 
condensation of the tightly-bound bosonic dimers made 
up of two fermions of opposite spin. In the crossover 
regime/unitary regime (l/feplosl ^ 1) bridging two well- 
understood regimes above, the gas is strongly interacting, 
which means that the interaction energy is comparable to 
the Fermi energy of the free Fermi gas of the same den- 
sity and defies the perturbative many-body techniques. 
Especially, the physics in the so-called "unitary limit" 
{\as\ = oo) is universal because only one remaining length 
scale, l/fcp, which is approximately the interparticle dis- 
tance, appears in the equation of state of the gas, while 
the interaction effect appears as a universal dimension- 
less parameter. Because the physics at unitarity is uni- 
versal, i.e., independent of its constituents and is non- 
perturbative due to the strong interaction, the unitary 
Fermi gas has been the test-bed for the various theoret- 
ical techniques developed up to now (see, e.g., Sec. VB 
and VC in Ref. [6]). 

Cold atomic gases also provide insights into other 
forms of matter. Inside neutron stars, there might be 
quark matter made of different ratios of quarks and this 
imbalance might result in new types of superfluid, such as 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase which 
has spatially non- uniform order parameter [10]. Even 
though the electric-charge neutrality of atoms might 
change the picture, we might get some ideas from the 
investigation on imbalanced atomic Fermi gases [11]. An- 
other example is the quark-gluon plasma (QGP) created 
in the Relativistic Heavy Ion CoUider (RHIC). QGP of 
temperature about 10^^ K was produced by smashing 
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nuclei together. Measurements on its expansion after its 
creation show that QGP is a nearly perfect fluid with 
very small shear viscosity. A cloud of atomic Fermi gases 
at unitarity shows the same strange behavior while the 
origins of the similarity between the hot and the cold 
matter are to be speculated [12]. 

Cold atoms in periodic potentials are more intrigu- 
ing because of their possible connection to the solid 
state/condensed matter physics (see, e.g., Refs. [5, 13] 
for reviews). The artificial periodic potential of light, 
so-called "optical lattice", is configured with pairs of 
counter-propagating laser fields which are detuned from 
atomic transition frequencies large enough that they act 
as practically frec-of-defect and conservative potentials 
that the atoms experience via optical dipole force. Geom- 
etry and depth of the lattice potential can be controlled 
by orienting the directions and changing the intensities 
of the laser beams, respectively. Moreover, macroscop- 
ically large lattice constant of the optical lattice com- 
pared with the lattice potential in solids is a great ad- 
vantage for experimental observation and manipulation; 
nowadays, in situ imaging and addressing of the strongly- 
correlated systems at the single-site level have been re- 
alized (e.g., Refs. [14]). The high controllability of both 
the optical lattice properties via laser fields and the in- 
teratomic interaction via the Feshbach resonance makes 
cold atomic gases serve as quantum simulators for vari- 
ous models of condensed matter physics. Many phenom- 
ena of solid sate/condensed matter physics have already 
been observed or realized in cold atomic gases on the op- 
tical lattices, including band structure, Bloch oscillation, 
Landau-Zcner tunneling, and superfluid-Mott insulator 
transition [15, 16]. The precision control and measure- 
ment at the level of one lattice site enable the cold atomic 
gases on optical lattice to be applied in inventing new ma- 
nipulation techniques and novel devices such as matter- 
wave interferometry, optical lattice clocks, and quantum 
registers [17, 18]. 

In this review article, we mainly explore the superfluid 
properties of cold atomic gases in the optical lattice at 
zero temperature. Superfluidity is the most well-known 
quantum phase of ultracold atomic gases and is prevalent 
in many other systems such as superconductors, super- 
fluid helium, and superfluid neutrons in "pasta" phases 
(see, e.g., Rcf. [19]) of neutron star crusts. Knowing its 
properties is also important both for judging whether the 
superfluid state is approached in experiments and apply- 
ing its properties in controlling the system. Depending 
on the physical conditions, cold atomic gases in optical 
lattices show various interesting phases other than the su- 
perfluid phase, for example, simulating Hubbard models 
and some other quantum phases will be explained very 
briefly in Sec. VI. In our discussion, we mainly consider 
atomic Fermi gases partly because weakly interacting 
atomic Bose gases can be considered as an extreme limit 
in the BCS-BEC crossover of the atomic Fermi gases and 
partly because physics of the cold Fermi gases is richer. 
Bose gases are discussed separately in the cases where we 



cannot find the corresponding phenomena in the Fermi 
gases. This does not mean that the cold Bose gases are 
less interesting and less important. The tools and tech- 
nique developed for creating the BEG were the building 
block for those for creating superfiuid Fermi gases and 
superfluid properties of the Bose gases might be more 
useful in some applications due to its higher superfluid 
transition temperature. 

In Sec. II, our system and the theoretical frameworks 
are presented. There, we explain the mean-field theory 
and the hydrodynamic theory depending on the validity 
regime of each theory. Using cold atomic gases, vari- 
ous parameter regimes become accessible and the mean- 
field theory in the continuum model is a powerful tool 
which allows us to study superfluid states covering such 
a wide region. Hydrodynamic theory, though its validity 
regime is limited, can provide precise analytic predictions 
in some limit, which can be served as complementary to 
the mean-flcld theory. Based on these methods, we study 
various aspects of the superfluid state in a lattice poten- 
tial periodic along one spatial direction. This setup is 
the simplest but contains essential effects of the lattice 
potential. 

Section HI discusses basic static and long-wavelength 
properties of the cold Fermi gases in the optical lattice. 
Speciflcally, the incompressibility and the effective mass 
are obtained from the equation of state (EOS). Focusing 
on the unitary Fermi g typical example, here we 

show that these properties in the lattice can be strongly 
modified from those in free space [29]. Formation of 
bound pairs is assisted by the periodic potential and this 
results in a qualitative change of the EOS. 

In Sec. IV, we examine the stability of the superfluid 
flow in the optical lattice. Mainly the energetic instabil- 
ity and the corresponding critical velocity of the super- 
fiow are discussed. Stability of the superflow is a fun- 
damental problem of the superfluidity and cold atomic 
gases allow us to study this important problem for vari- 
ous parameter regimes and for both the Bose and Fermi 
superfluids. Unlike an obstacle in uniform superfluids, 
periodic potential can modify the excitation spectrum 
and this manifests itself in the behavior of the critical 
velocity [27, 55]. 

In Sec. V, we investigate the energy band structures 
of cold atomic Fermi gases in an optical lattice along the 
BGS-BEG crossover. By tuning the interatomic inter- 
action strength, the nonlinear effect of the interaction 
energy can be made to dominate over external periodic 
potential. In this parameter regime, "swallowtail" en- 
ergy band structures emerge and the physical properties 
of the system is affected by their appearance [74]. This 
is one example that cold atomic gases can show physical 
phenomena unaccessible in other systems. 

In Sec. VI, we discuss quantum phases of cold atomic 
gases in optical lattices and show how they play their role 
as "quantum simulators" of Hubbard models. Before the 
advent of the cold atomic gases in optical lattices, the 
Hubbard model has been considered as an approximate 
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toy model of more complicated real system and been one 
of the central research topics in solid-state/condensed- 
matter physics. High-precision control of the cold atomic 
gases in optical lattices now makes the Hubbard models 
realized in experiments and we might address open ques- 
tions of quantum magnetism and high Tc superconduc- 
tivity. 

Finally, Sec. VH contains the summary and some per- 
spectives on the cold atomic gases in optical lattices. 



(LDA) is characterized by l/ryfino ^ {d/O"^ ^ 1 corre- 
sponding to the lattice with low fineness Tyane ^ 1- 

In the present article, we shall consider a large pa- 
rameter region covering weak to strong lattices, ryheight = 
0(10-2) - 0(10) (s - 0.1 - 5), and low to high fine 
lattices, rjunc 0.1 - 10. 



B. Mean- field theory in tlie continuum model 



II. THEORETICAL FRAMEWORKS 

A. Setup of the system and periodic potential 

In the present article, we discuss superfiuid fiow made 
of either fermionic or bosonic cold atomic gases subject 
to an optical lattice. For concreteness, we consider one of 
the most typical cases: the system is three dimensional 
(3D) and the potential is periodic in one dimension with 
the following form: 

Vcxt(r) = sEu sin^ qbz = Vq sin^ qbz. (1) 

Here, Vb = si?R is the lattice height, s is the lattice inten- 
sity in dimensionless units, Ej^ = h^q^/2m is the recoil 
energy, m is the mass of atoms, = tt/cI is the Bragg 
wave vector, and d is the lattice constant. For simplicity, 
we also assume that the supercurrent is in the z direction 
and thus the system is uniform in the transverse (i.e., x 
and y) directions. Throughout the present article, we set 
the temperature T = 0. 

Before giving a detailed description of the theoretical 
framework, it is useful to summarize the scales in this sys- 
tem. The periodic potential is characterized by two en- 
ergy scales: one is the recoil energy E^, which is directly 
related to the lattice constant d, and the other is the 
lattice height Vq = sEj^. Regarding Bose-Einstcin con- 
densates of bosonic atoms, a characteristic energy scale 
is the interaction energy gn, where g — 47rh?as/m is the 
interaction strength and n is the density of atoms. On 
the other hand, in the case of superfluid Fermi gases, 
the total energy is on the order of the Fermi energy 
Ey = h^kp/{2m) with the Fermi wave number k-p = 
(STT^n)^/"^ corresponding to the uniform non-interacting 
Fermi gas with the same density n. Therefore, the rel- 
ative effect of the lattice strength is given by the ra- 
tio 77hcight = Vo/gn for bosons and r?height = Vq/Ef for 
fcrmions. Likewise, relative fineness of the lattice (com- 
pared to the healing length) is characterized by the ratio 
?7fino ~ (s^Z-E'r)"^ (?/^)^ for bosons and, for fcrmions, 
7yfi„(, = (Ep/E^)^^ ~ (k^d)^'^, which is at most ~ {S,/dY 
near unitarity. Here, ^ is the healing length given as 
^ = h/{2mgny/^ for Bose superfluids and as ^ ~ kp^ 
for Fermi superfluids at unitarity, which is consistent 
with the BCS coherence length ^bcs = fivp/A, where 
vp = hkp/m and A is the pairing gap. We can also 
say that the validity of the local density approximation 



To study the superfluidity in such various regimes in a 
unified manner, one of the most useful theoretical frame- 
works is the mean-field theory in the continuum model. 
Since our system is a superfluid in three dimensions and 
the number of particles in each site is inflnite in our setup, 
effects of quantum fluctuations, which are not captured 
by the mean-field theory, may be small. It is also noted 
that widely-used tight-binding models are invalid in the 
weak lattice regime of ?7hcight ^ 1- 

For dilute BECs at zero temperature, the system is 
well described by the Gross-Pitacvskii (GP) equation [4, 
20, 21]: 

-|^5f* + yext(^)* + 5l*P*-M*, (2) 
Zra 

where ^'(z) = \/ n{z) exp[i0(z)] is the condensate wave 
function, is its phase, and fj, is the chemical po- 

tential. The local superfiuid velocity is given by v{z) = 
{h/m)dz(j>{z). 

For superfiuid Fermi gases, we consider a balanced sys- 
tem of attractively interacting (pseudo)spin 1/2 fcrmions, 
where the density of each spin component is n/2. To de- 
scribe the BCS-BEC crossover at zero temperature, we 
use the Bogohubov-de Gennes (BdG) equations [6, 22], 

fH'ir) A(r) \ f u,{r) \ _ / u,(r) \ 
l,A*(r) -H'{r) )[u,{r) )^'^[v,{r) ) ' 

with H'{y) — —h^V^ /2m + Vc^t ~ Mi(r) and Wi(r) are 
quasiparticle wave functions, and Ci is the correspond- 
ing quasiparticle energy. The pairing field A(r) and the 
chemical potential ^ in Eq. (3) are self-consistently de- 
termined from the gap equation 

^{v) = -gY,u,{v)v*{v) , (4) 

i 

together with the constraint on the average number den- 
sity 

^=|E/ KWl'rfr^l jn{v)dv, (5) 

with n(r) = 2 ^ - |t;.j(r)p. Here g is the coupling constant 
for the contact interaction and V is the volume of the 
system. The average energy density e can be calculated 
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as 



e=- Idv 



2m 



9 



|A(r)|^ 



(6) 



For contact interactions, the right hand side of Eq. (4) 
has uhraviolet divergence, which has to be regularized by 
replacing the bare coupling constant g with the two-body 
T- matrix related to the s-wave scattering length [23-26] . 
A standard scheme [23] is to introduce the cutoff energy 
Ec = /i^fcp/2m in the sum over the BdG eigenstates and 
replace g by the following relation: 



1 

9 

(0) = fi2L2, 



E 

k<kc 



1 



2e 



(7) 



with e^''^ = h'e/2m. 

In the presence of a supercurrent with wavevector Q = 
P/h (P is the quasimomentum for atoms rather than for 
pairs, thus it is defined in the range of |P| < hqB /2) in the 
z direction, one can write the quasiparticle wavefunctions 
in the Bloch form as Ui{r) ~ Wi(z)e*'^^e*'' '' and Vi{r) 
Vi{z)e' 



-iQz 



leading to the pairing field as 
A(r) = e*2'3^A(z). 



(8) 



Here A(z), Ui{z), and Vi{z) are complex functions with 
period d and the wave vector i\kz\ < ^b) lies in the 
first Brillouin zone. This Bloch decomposition trans- 
forms Eq. (3) into the following BdG equations for Uiiz) 
and Vi{z) : 



Hq{z) K{z) 
A*(z) -H^q{z) 



Ui{z) 



Ui{z) 



(9) 



where 



HQiz) 



t2 r 



2m 



kl + {-id, + Q + k,Y 



Here, k\ = k'^ + ky and the label i represents the wave 
vector k as well as the band index. 



C. Hydrodynamic theory 

When the local density approximation (LDA) is valid, 
such that the typical length scale of the density varia- 
tion given by d is much larger than the healing length ^ 
of the superfluid, hydrodynamic theory in the LDA can 
be useful [27]. In hydrodynamic theory, we describe the 
system in terms of the density field n{z) and the (quasi- 
)momentum field P{z) [or the velocity field v{z)]. The 
LDA assumes that, locally, the system behaves like a uni- 
form gas; thus the energy density e(n, P) can be written 
in the form 



e(n, P) = nP^/2m + e(n, 0) 



(10) 



and one can define the local chemical potential /x(n) = 
de{n,0)/dn. The density profile of the gas at rest in 
the presence of the external potential can be obtained 
from the Thomas-Fermi relation: /.io = ^[n(z)] -I- Vcxtiz). 
If the gas is flowing with a constant current density 
j — n{z)v{z), the Bernoulli equation for the stationary 
velocity field v{z) is 



] 



n{z) 



(11) 



where is the z-independcnt value of the chemical po- 
tential. 

In the following two typical cases, the uniform gas has 
a polytropic equation of state: 



(12) 



1) A dilute Bose gas with repulsive interaction, where 
7 = 1 and a = g = ATrfi^as/m.. 2) A dilute Fermi 
gas at unitarity, where fj,{n) = (1 + (3)Ep — [(1 + 
/3)(37r2)2/3s2/(2m)]n2/3, i.e., 7 = 2/3 and a = (1 + 
/3)(37r2)2/'^fi,2/(2?7i). Here /3 is a universal parameter, 
which is negative and its absolute value is of order unity 
accounting for attractive interatomic interactions [6, 28]. 
Using the equation of state one can write 



mcliz) = n^n{n) = -fii{n) , 



(13) 



where Cs{z) is the local sound velocity, which depends 
on z through the density profile n{z). In a uniform 

gas of density n, the sound velocity is given by cf' = 
[7^(n)/m]^/2. 



III. EQUATION OF STATE, 
INCOMPRESSIBILITY, AND EFFECTIVE MASS 

In this section, focusing on superfluid Fermi gases at 
unitarity, we discuss effects of the periodic potential on 
macroscopic and static properties of the fluid such as 
equation of state, incompressibility, and effective mass 
[29]. The important point is that the periodic poten- 
tial favors the formation of bound molecules in a two- 
component Fermi gas even at unitarity [30] (see also, 
e.g., Refs. [31, 32]). The emergence of the lattice-induced 
bound states drastically changes the above macroscopic 
and static properties from those of uniform systems in 
the strong lattice regime of ryhoight ^ 1 [29]. Such an 
effect is absent in ideal Fermi gases and BECs of repul- 
sively interacting bosonic atoms which can be considered 
as two limits in the BCS-BEC crossover. 



A. Basic equations 

At zero temperature, the chemical potential fi (the 
equation of state) and the current j of a superfluid Fermi 
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gas in a lattice are given by the derivative of the aver- 
age energy density e = E/V with respect to the average 
(coarse-grained) density n [33] and the average quasimo- 
mentum P of the bulk superflow, respectively (hereafter, 
for notational simplicity, we omit " ~ " for the coarse- 
grained quantities, which should not be confused with 
the local quantities): 



/i 



de{n,P) 



J 



de{n,P) 



(14) 



dn ' dP 

The incomprcssibility (or inverse compressibility) 
and the effective mass m* , are given by the second deriva- 
tives of e with respect to n and P: 



K ^ =n 



d'^ein, P) dnin, P) 
~ n 



1 ld'^e{n,P) 



dn ' 
iaj(n,P) 



(15) 
(16) 



m* n dP^ n dP 
We calculate these quantities for P = 0, i.e. for a gas at 
rest, in the periodic potential. 

In the absence of the lattice potential (s = 0), thermo- 
dynamic properties of unitary Fermi gases show a univer- 
sal behavior: the only relevant length scale is the inter- 
particlc distance fixed by k-p . Due to translational invari- 
ance, one can write e{n,P) = e(n,0) -I- nP'^/2m so that 
j = nP/m and m* = m. Furthermore, the energy density 
at P = can be written as e(n, 0) = (l-|-/3)e°(n, 0), where 
e°(n, 0) = (3/5)n£'F oc n^/"^ is the energy density of the 
ideal Fermi gas. Thus we have yU = (1 -f /3)i?F + P^ /2m 
and K'^ = {2/2,){l + l3)Ep. 



B. Equation of state and density profile 

When the lattice height s is large, the periodic po- 
tential favors the formation of bound molecules. In the 
strong lattice limit ryhoight = sE^/E-p ^ 1, the system 
tends to be a BEC of lattice-induced bosonic molecules. 
Therefore, in this regime, the chemical potential shows 
a linear density dependence, fj. oc n, as shown by the 
red solid line in the inset of Fig. 1 calculated for uni- 
tary Fermi gases in a lattice with s — 5. This is clearly 
different from the density dependence of the chemical po- 
tential in the uniform system (s = 0), ^ oc n^/^, as shown 
by the blue dashed line in the same inset. It is also noted 
that, for s » 1, this linear density dependence persists 
even at relatively large densities so that Ep/sEj^^ ^ 1 
(e.g., nq^^ = 0.1 corresponds to Ey/E-r ~ 2.1 and thus 
E-p/sEn ~ 0.41 in the case of this figure) because the 
system effectively behaves like a 2D system in this den- 
sity region due to the bandgap in the the longitudinal 
degree of freedom. 

This drastic change of EOS manifests itself as a change 
of the coarse-grained density profile when a harmonic 
confinement potential Vho(r) is added to the periodic po- 
tential. The coarse-grained density profile, n(r), is cal- 
culated using the LDA for /.i: 
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fio = M["(r)] + Vho(r)- 



(17) 



FIG. 1: (Color online) Coarse-grained density profiles of a 
trapped unitary Fermi gas, n{r± = 0, z) for s = and 5 in 
units of the central density n(0) — 0.0869^1 calculated for 
s — (this local density corresponds to Ep/En — 1.88). The 
quantity i?!"' is the axial Thomas- Fermi radius for s — 0. The 
inset shows the density dependence of the chemical potential 
of unitary Fermi gases. Here fio is the chemical potential in 
the limit of n = 0. This figure is taken from Ref. [29] 



Here fi[n\ is the local chemical potential as a function of 
the coarse-grained density n obtained by the BdG cal- 
culation for the lattice system and /io is the chemical 
potential of the system fixed by the normalization con- 
dition J drn{r) = N. Figure 1 clearly shows that, for 
s = 5, the profile takes the form of an inverted parabola, 
reflecting the linear density dependence of the chemical 
potential (see inset). In this calculation, we assume an 
isotropic harmonic potential, Vho(r) = muj'^r'^/2, where 
Lo is the trapping frequency, fnu/En = 0.01, and the num- 
ber of particles N = 10^; these parameters are close to 
the experimental ones in Ref. [34]. 



C. Incompressibility and effective mass 

The formation of molecules induced by the lattice has 
important consequences also on and m*. Due to 
the linear density dependence of the chemical potential 
in the strong lattice regime (jyhcight ^ 1 or E-p/E^i ^ s; 
or low density limit for a fixed value of s), is also 
proportional to n and k~^/k~^(s = 0) oc n^/"^ — > for 
Ef/E^^ — ^ [sec Fig. 2(a)]. This means that the gas 
becomes highly compressible in the presence of a strong 
lattice. This is in a strong contrast to the ideal Fermi 
gas corresponding the the BCS limit, which gives nonzero 
values of k^^ / k^^{s = 0) ~ 1 even in the same limit (see 
Fig. 2 in Ref. [29]). On the other hand, in the weak lat- 
tice limit (or high density limit for a fixed value of s), 
the system reduces to a uniform gas. By using an hydro- 
dynamic theory, which is valid when Ep / E-^ ^ 1 , and 
expanding with respect to the small parameter sE-n/E-p, 
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we obtain k 

[291 



of unitary Fermi gases in this regime as 



1 , 



P) 



Ef 

0[{sEk/E_ 



(18) 



This is shown by dotted lines in Fig. 2(a). Note that, in 
this regime, / k~^{s = 0) > 1 and it decreases to unity 
with increasing sEji/Ep. Therefore, k^^/k^^{s = 0) 
should take a maximum value larger than unity in the 
intermediate regime of E^/E-n ~ 1 as can be seen in 
Fig. 2(a), which is mainly caused by the bandgap in the 
longitudinal degree of freedom. 

Since the tunneling rate between neighboring sites, 
which is related to the (inverse) effective mass is 
exponentially suppressed with the mass, the formation of 
molecules induced by the lattice can yield a drastic en- 
hancement of TO* for s ^ 1 in the low density limit [Fig. 
2(b)]. This enhancement makes to* much larger com- 
pared to the ideal Fermi gases (see Fig. 2 in Ref. [29]) 
and to the BECs of repulsively interacting Bose gases 
(see Fig. 4 in Ref. [35]) with the same mass to. As 
Ep / E^ increases, the effective mass exhibits a maximum 
at Ep / En ~ 1 due to the bandgap in the transverse 
degree of freedom, then it decreases to the bare mass, 
771* = TO. Hydrodynamic theory can explain the behav- 
ior of TO* of unitary Fermi gases for small sEn/E-p [29]: 



TO 
TO 



9 

32 



sE^ 
Ep 



^0[{sEp,/Epf] . (19) 



The numerical factor in the second term shows that the 
effect of the lattice is stronger for to* than for k^^. It 
is worth comparing the results with the case of bosonic 
atoms, where to* decreases monotonically with increas- 
ing density since the interaction broadens the condensate 
wave function and favors the tunneling [35] . 

In Fig. 3, we show the sound velocity of the unitary 
Fermi gases in a lattice. 



K 

TO* 



(20) 



calculated from the above results of and to*. It shows 
a significant reduction compared to the uniform system 
mainly due to the larger effective mass to* /to > 1 except 
for the low density (more precisely, strong lattice) limit in 
which and thus Cg show an abrupt reduction. Using 
Eqs. (18) and (19), we obtain the expression of the sound 
velocity in the weak lattice limit as 



+ +0[{sEn/Epr 



(21) 

where ci°^ = [{2/3){l + (3)Ep /m]^^'^ is the sound velocity 
for the uniform system. 
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FIG. 2: (Color online) Incompressibility and eflective 
mass m* of unitary Fermi gases for s = 1 (red), 2.5 (blue), 
and 5 (green). Asymptotic expressions (18) and (19) obtained 
by the hydrodynamic theory are shown by the dotted lines. 
Open circles in panel (b) show m* obtained in Ref. [30] by 
solving the Schrodinger equation for two-body problem. The 
s = 1 results for m* are also shown in the inset in the linear 
scale. This figure is adapted from Ref. [29]. 
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FIG. 3: (Color online) Sound velocity Cs of unitary Fermi 
gases in a lattice in units of the sound velocity ci'^^ — 
[(2/3)(l + /3)EF/mY/'^ for the uniform system. As in Fig. 
2, red, blue and green lines correspond to s = 1, 2.5, and 5, 
respectively. This figure is taken from Ref. [29] 
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IV. STABILITY 
A. Landau criterion 

Stability of superfluid flow is one of the most funda- 
mental issues of superfluidity, which was pioneered by 
Landau [36] . He predicted that there is a critical veloc- 
ity Vc of the superflow above which the kinetic energy of 
the superfluid is large enough to be dissipated by creat- 
ing excitations (see, e.g., Refs. [36-39]). This instability 
is called Landau or energetic instability and its critical 
velocity is Landau critical velocity. 

The celebrated Landau criterion for energetic instabil- 
ity of uniform superfluids is given by [36-39] 

v>v^ = min i^-y^ , (22) 

where v is the velocity of the superflow, e(p) is the exci- 
tation spectrum in the static case (u = 0), and p is the 
magnitude of momentum p of an excitation in the co- 
moving frame of the fluid. Here, is determined by the 
condition for which there starts to exist a momentum p 
at which the excitation spectrum in the comoving frame 
of a perturber (it can be an obstacle moving in the fluid 
or a vessel in which the fluid flows) is zero or negative. 

For superfluids of weakly interacting Bose gases, ex- 
citation spectrum is given by the Bogoliubov dispersion 
relation (e.g., Refs. [37-40]), 

(23) 

with 

cs ^ .[^ (24) 

being the sound velocity defined by Eq. (20) [note that 
ji = gn and thus = gn, and m* = m for uniform 
BECs]. Thus, from Eqs. (22) and (24), it is easily seen 
that the Landau critical velocity is given by the sound 
velocity Vc = Cg forp = 0. This means that, in superfluids 
of dilute Bose gases, the energetic instability is caused by 
excitations of long-wavelength phonons. It is also noted 
that BECs of non-interacting Bose gases cannot show 
the superfluidity in a sense that Wc = and they cannot 
support a superflow. 

In superfluid Fermi gases, there is another mecha- 
nism which can cause the energetic instability: fermionic 
pair-breaking excitations. In the mean-field BCS theory, 
quasiparticle spectrum of uniform superfluid Fermi gases 
is given by 

e(.) = . (25) 



Thus the Landau critical velocity due to the pair- 
breaking excitations is given by [41] 

V, - ^^-(v/^I^TA^-Ai) . (26) 

In the deep BCS regime, where /i ~ E-p ^ A, we obtain 
Vc — A/pF with pf = hkp. 

In the BCS-BEC crossover of superfluid Fermi gases, 
where the above two kinds of excitations exist, the Lan- 
dau critical velocity is determined by either of them 
which gives smaller Vc- In the weakly interacting BCS 
regime (l/fepOs < —1), the pairing gap is exponen- 
tially smah, A ^ i^Fe^'^/^'^^l^^l, and Vc is set by the 
pair-breaking excitations. On the other hand, in the 
BEC regime (l/fepfls ^ 1): where the system consists of 
weakly interacting bosonic molecules, creation of long- 
wavelength superfluid phonon excitations causes the en- 
ergetic instability. In the unitary regime both mecha- 
nisms are suppressed and the critical velocity shows a 
maximum value [41-43]. 



B. Stability of superflow in lattice systems 

1. Energetic instability and dynamical instability 

For the onset of the energetic instability, energy dis- 
sipation is necessary in general; i.e., in closed systems, 

V > Vc is not a sufficient condition for a breakdown of 
the superflow and the flow could still persist even at 

V > Vc- The energetic instability corresponds to the sit- 
uation in which the system is located at a saddle point 
of the energy landscape (i.e., there is at least one direc- 
tion in which the curvature of the energy landscape is 
negative) . 

In the presence of a periodic potential, another type 
of instability called dynamical (or modulational) insta- 
bility can occur in addition to the energetic instability. 
The dynamical instability means that small perturba- 
tions on a stationary state grow exponentially in the pro- 
cess of (unitary) time evolution without dissipation. Sim- 
ilar to energetically unstable states, dynamically unsta- 
ble states are also located at saddle points in the energy 
landscape, but a difference from energetically unstable 
states is that there are kinematically allowed excitation 
processes which satisfy the energy and (quasi)momentum 
conservation. This means that energetic instability is a 
necessary condition for dynamical instability (see, e.g., 
Refs. [44, 45] for bosons and [46] for fermions), and there- 
fore, critical value of the (quasi)momentum for the dy- 
namical instability should be always larger than (or equal 
to) that for the energetic instability. Due to such a rea- 
son, we shall focus on the energetic instability hereafter 
[47]. 
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2. Determination of the critical velocity 

If the energetic instability is caused by long- wavelength 
superfluid phonon excitations, the critical velocity can 
be determined by hydrodynamic analysis of the excita- 
tions [27, 38, 52-55] (this should not be confused with 
the LDA hydrodynamics discussed in Sec. II C). This 
analysis is valid provided that the wavelength of the ex- 
citations which trigger the instability is much larger than 
the typical length scale of the density variation, i.e., the 
lattice constant d. 

We consider the continuity equation and the Euler 
equation for the coarse-grained density n and the coarse- 
grained (quasi)momentum P averaged over the length 
scale larger than the lattice constant. 



dn „ . dn d de 
h V • 7 = 1 = 

dt dt dzdP 



dP 

'dt 



dP d de 



dt dz dn 



= 0, 



(27) 
(28) 



where e = e(n, P) is the energy density of the super- 
fluid in the periodic potential for averaged density n 
and averaged (quasi)momentum P. Linearizing with 
respect to the perturbations of n{z,t) = uq + dn{z,t) 
and P{z,t) = Po+ dP{z,t) with Sn{z,t) and 5P{z,t) oc 
^iqz-iuit^ we obtain the dispersion relation of the long- 
wavelength phonon. 



d\ 
dndP 



q 



d'^e d'^e 
d^dP^ 



(29) 



Here, Huj and q arc the energy and the wavenumbcr of 
the excitation. In the first term (so-called the Doppler 
term), dndpe > 0. Thus the energetic instability occurs 
when uj{q) for (jf = — jgl becomes zero or negative: 



d^e 
dndP 



> 



d'^e d'^e 
dn? dP-^ 



(30) 



Using this condition, we determine the critical quasimo- 
mcntum Pc at which the equality of Eq. (30) holds, and 
finally we obtain the critical velocity from 



de 
'dP 



(31) 



It is noted that, for calculating Pc and Vc using 
Eqs. (30) and (31), what we only need is the energy 
density of the stationary states as a function of n and 
P. This can be obtained by solving, e.g., GP or BdG 
equations for the periodic potential. 

If the energetic instability of Fermi supcrfluids is 
caused by the pair-breaking excitations, the critical ve- 
locity can be determined by the quasiparticle energy 
spectrum obtained from the BdG equations. The en- 
ergetic instability by the pair-breaking excitations occurs 
when some quasiparticle energy becomes zero or neg- 
ative: 



We obtain a critical velocity for the pair-breaking exci- 
tations from Eq. (31) evaluated at the critical quasimo- 
mcntum determined by this condition. 



3. Critical velocity of superfluid Bose gases and superfluid 
unitary Fermi gases in a lattice 

First, we consider the situation where the LDA is valid, 
i.e., the lattice constant d is much larger than the heal- 
ing length ^. This condition corresponds to gn/Eji ^ 1 
for superfluid Bose gases and Ep/E^ ^ 1 for superfluid 
Fermi gases at unitary (see discussion in Sec. II A). In 
the framework of the LDA hydrodynamics explained in 
Sec. lie, the system is considered to become energeti- 
cally unstable if there exists some point z at which the 
local superfluid velocity v{z) is equal to or larger than 
the local sound velocity Cs{z). Suppose the external po- 
tential has a maximum at z = zq [i.e., V{zo) = Vq for our 
periodic potential], then at the same point the density 
is minimum, Cs{z) is minimum, and v{z) is maximum 
due to the current conservation, i.e., j — n(z)v{z) be- 
ing constant. This means that the superfluid becomes 
flrst unstable at z = zq. Using Eq. (13), the condition 
for the occurrence of the instability can be written as 
m[jc/nc{zo)]'^ = 7M[nc(zo)] = janj^zo), where n^z) is 
the density profile calculated at the critical current [56]. 
By inserting this condition into Eq. (11), one can obtain 
the following implicit relation for the critical current [27]: 



Jc 



1 



ma 



2/7 



2^ij 



■7 



^fo) 



1 4+1 



(33) 



< 0. 



(32) 



It is worth noticing that this equation contains only z- 
independent quantities. It is also independent of the 
shape of the external potential: the only relevant pa- 
rameter is its maximum value V{zq). Moreover, it can 
be applied to both Bose gases and unitary Fermi gases 
(see also Refs. [57-59] for bosons). 

In Fig. 4, we plot the critical velocity in a lattice ob- 
tained from the hydrodynamic expression (33) for BECs 
[panel (a)] and unitary Fermi gases [panel (b)] by the 
thick solid red lines. In both cases, the critical velocity 
f c = jc/no (no is the average density) is normalized to the 
value of the sound velocity in the uniform gas, c"s \ and is 
plotted as a function of Vo/Mj=o- The limit Vq/ ^j=o ~^ 
corresponds to the usual Landau criterion for a uniform 
supcrfiuid fiow in the presence of a small external pertur- 
bation, i.e., a critical velocity equal to the sound velocity 
of the gas. In this hydrodynamic scheme, as mentioned 
before, the critical velocity decreases when Vq increases 
mainly because the density has a local depletion and the 
velocity has a corresponding local maximum, so that the 
Landau instability occurs earlier. When Vq ~ /ij^, the 
density exactly vanishes at z = zq and hence the critical 
velocity goes to zero. 

Next, we discuss the critical velocity of the energetic 
instability beyond the LDA [27] . Here we use the energy 
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FIG. 4: (Color online) Critical velocity Vc for energetic in- 
stability of superfluids in a ID periodic potential. Panel (a) 
is for superfluids of dilute Bose gases and panel (b) for super- 
fluids of dilute Fermi gases at unitarity. The critical velocity 
is given in units of the sound velocity of a uniform gas, ci'^\ 
and is plotted as a function of the maximum of the exter- 
nal potential in units of the chemical potential fJ,j=o of the 
superfluid at rest. Thick solid lines: prediction of the hydro- 
dynamic theory within the LDA, calculated from Eq. (33). 
Symbols: results obtained from the numerical solutions of 
the GP equation [panel (a)] and the BdG equations [panel 
(b)]. The thinner black solid lines are the tight-binding pre- 
diction (35). Dashed lines are guides to the eye. This figure 
is adapted from Ref. [27]. 



density e(n,P) of superfluids in a periodic potential cal- 
culated by the mean-field theory in the continuum model, 
i.e., the GP equation for bosons and the BdG equations 
for fermions. Based on the energy density, we determine 
the critical quasimomcntum and the critical velocity from 
Eqs. (30) and (31). 

For Bose superfluids, we plot Vc for various values of 
gn/EB, in Fig. 4(a). We can clearly see that these re- 
sults approach the LDA prediction for gn/E^^ ^ 1, as 
expected. We also note that Vc exhibits a plateau for 
gn/Eii < 1 (i.e., £, ^ d) and small Vq. This can be 
understood as follows. If the healing length ^ is larger 



than the lattice spacing d and Vo/gri is not too large, the 
energy associated with quantum pressure, which is pro- 
portional to 1/^^, acts against local deformations of the 
order parameter and the latter remains almost unaffected 
by the modulation of the external potential. This is the 
regime of the plateau in Fig. 4. In terms of Eq. (30), this 
regime occurs when the left-hand side is ~ P/m and the 

right-hand side is ~ c's^ , so that the critical quasimomcn- 
tum obeys the relation Pc/m ~ cf\ which is the usual 
Landau criterion for a uniform superfluid in the presence 
of small pcrturbcrs. By increasing Vb, this regime ends 
when /i ~ E'r. 

If we further increase Vb , the chemical potential /i be- 
comes larger than En, the density is forced to oscillate, 
and Vc/cf^ starts to decrease. The system eventually 
reaches a regime of weakly coupled superfluids separated 
by strong barriers, which is well described by the tight- 
binding approximation (also for gn/E^ > 1, the system 
enters this regime when Vq is sufficiently large). There, 
the energy density is given by the sinusoidal form with 
respect to P as 

e{n, P) ^ e{n, 0) + Sj [1 - cos (^P/Pedgc)] • (34) 

Here, Podge is the quasimomcntum at the edge of the 
first Brillouin zone, i.e., Pedge = ^^b for superfluids of 
bosonic atoms (Podge = ^9b/2 for those of fermionic 
atoms). The quantity Sj ~ nP^^^^/ir'^rn* corresponds to 
the half width of the lowest Bloch band. Because of the 
sinusoidal shape of the energy density and large effective 
mass in the tight-binding limit, the critical quasimomen- 
tuni is around Podge/2. Thus we sec that, from Eq. (31), 
the critical velocity is determined by the effective mass 
as 

vc^-^. (35) 
TT m 

The values of Vc obtained from Eq. (35), with m* ex- 
tracted from the GP calculation of e{n,P), are plotted 
in Fig. 4(a) by thinner black solid lines for gn/Ej^ = 0.4 
and 1 in the region of Vo/nj=o }t 2. 

Wc calculate the critical velocity also by another 
method based on the complete linear stability analysis 
for the GP energy functional as in Refs. [44, 52, 60] (see 
also Refs. [38, 45]). We have checked that the results 
agree with those obtained from Eq. (30) based on the 
hydrodynamic analysis within 1% in the whole range of 
gn/E^ and Vb considered in the present work. This con- 
firms that the energetic instability in the periodic poten- 
tial is triggered by long-wavelength excitations, since the 
excitation energy of the sound mode is the smallest in 
this limit. 

For superfluid unitary Fermi gases, we plot Vc for var- 
ious values of Ep/Eji in Fig. 4(b). We observe qual- 
itatively similar results to those for bosons plotted in 
Fig. 4(a). For E-p/En < 1, the critical velocity Vc shows 

a plateau at small Vb and it decreases from ~ ci"^ with 
increasing Vq. For larger E'f/Pr > 1, Vc approaches the 
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FIG. 5: (Color online) Critical velocity Vc of the energetic 
instability for Ep/Eji = 0.1 with s = 1 and 5 in the BCS- 
BEC crossover. Open circles and filled squares show the crit- 
ical velocity due to long-wavelength phonons and fermionic 
pair-breaking excitations, respectively. The horizontal dotted 
line represents the value of the sound velocity ci°' of uniform 
system at unitarity, ci°V'fF = (1 + P^^^/VS ~ 0.443. The 
red solid line is a guide to the eye. 



LDA result (thick solid red line) . In the regime of large 
Vq such that Vq/a* is sufficiently large, Vc is well described 
by the tight-binding expression (35) plotted by thinner 
black solid lines. Here we use m* calculated from BdG 
equations. We also note that the pair-breaking excita- 
tions are irrelevant to the energetic instability at unitar- 
ity except for very low densities such that Ep/Ej^^ <C 1 
and small but nonzero Vq. These excitations are impor- 
tant in the BCS side of the BCS-BEC crossover discussed 
in See. IV B 4. 

Finally, we would like to point out that, in the LDA 
limit, while the quantities discussed in Sec. Ill approach 
the value in the uniform system, the critical velocity of 
the energetic instability is most strongly affected by the 
lattice. Figure 4 clearly shows that, as we increase the 
lattice height, the critical velocity decreases from that in 
the uniform system most rapidly in the LDA limit. The 
main reason of this opposite tendency is that the critical 
velocity is determined only around the potential maxima 
while the quantities discussed in Sec. Ill are determined 
by the contribution from whole region of the system. 



4- Along the BCS-BEC crossover 

Here we extend our discussion on the Landau critical 
velocity to the BCS-BEC crossover region [55]. As in 
the uniform systems, both the long-wavelength phonon 
excitations and pair-breaking excitations can be relevant 
to the energetic instability depending on the interaction 
parameter l/fepfls- However, there is an additional effect 
due to the lattice: when the lattice height is much larger 
than the Fermi energy, the periodic potential can make 
pairs of atoms to be strongly bound even in the BCS 
regime, and the pair-breaking excitations are suppressed 
[55]. 

In Fig. 5, we plot Vc for superfluid Fermi gases in 
the BCS-BEC crossover for different values of the lattice 
height with s = 1 and 5. Here we set Ep/En = 0.1 as an 
example. The open circles show Vc of the energetic insta- 
bility caused by the long-wavelength superfluid phonon 
excitations and the filled squares show that by the pair- 
breaking excitations. For moderate lattice strength of 
s = 1, the result of Vc is qualitatively the same as that of 
the uniform system. In the BCS side, smallest Vc is given 
by the pair-breaking excitations while, in the BEC side, 
Vc is set by the long-wavelength phonon excitations, and 
around unitarity Vc shows a maximum. For a larger value 
of .s = 5, however, the result is very different from the 
uniform case. Due to the lattice-induced binding, pair 
breaking is suppressed so that it does not cause the ener- 
getic instability even in the BCS side (Note that there is 
no filled square for s — 5. This means that we do not have 
negative quasiparticle energy for any value of P in the 
whole Brillouin zone. This point will be discussed later.). 
Throughout the calculated region of — 1 < l/k-pUs < 1, 
Vc is determined only by the long- wavelength phonon ex- 
citations and Vc decreases monotonically with increasing 
l/kpUs- 

The effect of the lattice-induced molecular formation 
can be clearly seen in the quasiparticle energy spectrum 
ei and in the enhancement of the order parameter A(2). 
In Fig. 6, we show and |A(z)| in the case of s = 5, 
E-p/E-R, = 0.1, and \/kpas = —1. The spectrum for 
P = shows a quadratic dependence of with a pos- 
itive curvature around fc^ = and there are no minima 
at fcj ^ 0. Even though the figure represents a case in 
the BCS regime, the structure of is consistent with 
the formation of bound pairs. We can also see that the 
spectrum never becomes negative for any values of P in 
the first Brillouin zone. In the inset of the same figure, 
we show the amplitude | A(z)| of the order parameter at 
P = 0. Here we see a large enhancement of |A(z)|, near 
z = 0, compared to the uniform system, which shows 
the formation of bosonic bound molecules. We also note 
that the minimum value of |A(z)| at z/d = ±1 is smaller 
than, but still comparable to the value of |A| in the uni- 
form case, suggesting that the system is indeed in the 
superfluid phase. 
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FIG. 6: (Color online) Lowest band of the quasiparticle en- 
ergy spectrum for large lattice height with s — 5 and 
Ep/En = 0.1 (i.e., Ef/Vq = 0.02) in the BCS regime at 
l/kpas — —1. Here, we show the first radial branch with 
k\ = + ky = 0, which always gives the smallest values of 
£; in this case. The inset shows the amplitude |A(z)| of the 
order parameter at P = 0. The horizontal dotted line shows 
the amplitude of the order parameter for the uniform system 
at the same values of l/kpa^ = —1. This figure is adapted 
from Ref. [55]. 



ity was first experimentally studied in Ref. [66]. In this 
experiment, they used a cigar-shaped BEC which under- 
goes the center-of-mass oscillation in a harmonic trap in 
the presence of a weak ID optical lattice, and measured 
the critical velocity. Although their original conclusion 
was that they have measured the Landau critical velocity 
for the energetic instability, further careful experimental 
[67] and theoretical [44, 60] follow-up studies have clar- 
ified that the instability observed in Ref. [66] was a dy- 
namical instability. In this follow-up experiment by De 
Sarlo et al. [67], they employed an improved setup, i.e., 
a ID optical lattice moving at constant and tunable ve- 
locities instead of an oscillating BEC in a static lattice. 
With this new setup, they succeeded in observing both 
energetic [67] and dynamical [67, 68] instabilities, and ob- 
tained a very good agreement with theoretical predictions 
[60, 67, 68]. It is worthwhile to stress that this under- 
standing was finally obtained by a continuous effort over 
a few years by the same group. 

Experimental study on stability of superfluid Fermi 
gases in optical lattices was carried out by Miller et al. 
[34]. Similar to the experiment of Ref. [67], they also 
used a ID lattice moving at constant and tunable ve- 
locities. A different point is that, instead of imposing a 
periodic potential on the whole cloud, they produced a 
lattice potential only in the central region of the cloud. 
They measured the critical velocity of the energetic in- 
stability in the BCS-BEC crossover and found that it 
is largest around unitarity. They also did a systemati- 
cal measurement of the critical velocity at unitarity for 
various lattice height. However, there is a significant dis- 
crepancy from a theoretical prediction [27] and further 
studies are needed. 



5. Experiments 

In closing this section, we briefly summarize experi- 
mental studies on stability of the superfluid. Using cold 
atomic gases, stability of the superfluid and its critical 
velocity was first experimentally studied in Ref. [61] and 
further examined in Ref. [62]. In these experiments, they 
put a large (diameter ^ ^) and strong (height ^ fi) vi- 
brating circular potential in a BEC. However, it has been 
concluded that what they observed in these works were 
not likely to be the energetic instability: dynamical nu- 
cleation of vortices by vibrating potential [62] . Recently, 
Ramanathan et al. performed a new experiment on sta- 
bility of the superflow with a different setup [63]. They 
used a BEC flowing in a toroidal trap with a tunable 
weak link (width of the barrier along the flow direction is 
much larger than ^). In this experiment, they obtained 
the critical velocity consistent with the energetic insta- 
bility due to the vortex excitations [64]. For 2D Bose 
gases, their stability of the superfluidity and the critical 
velocity was also studied recently [65] . 

Regarding the superflow in optical lattices, its stabil- 



V. ENERGY BAND STRUCTURE 

In this section we discuss how the superfluidity can af- 
fect the energy band structures of ultracold atomic gases 
in periodic optical potentials. Starting with the nonin- 
teracting particles in a periodic potential, we obtain the 
well-known sinusoidal energy band structure. When the 
interparticle interaction is turned on in such a way that 
the superfiuidity appears, we see a drastic change in en- 
ergy band structure. For BECs, it has been pointed out 
that the interaction can change the Bloch band struc- 
ture drastically, causing the appearance of a loop struc- 
ture called "swallowtail" in the energy dispersion [69-72] . 
This is due to the competition between the external pe- 
riodic potential and the nonlinear mean-field interaction: 
the former favors a sinusoidal band structure, while the 
latter tends to make the density smoother and the energy 
dispersion quadratic. When the nonlinearity wins, the ef- 
fect of the external potential is screened and a swallowtail 
energy loop appears [73]. This nonlinear effect requires 
the existence of an order parameter and, consequently, 
the emergence of swallowtails can be viewed as a pecu- 
liar manifestation of superfluidity in periodic potentials. 
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FIG. 7: (Color online) (a) Energy E per particle as a function 
of the quasimomentum P for various values of 1/fcpas; (b) 
half-width of the swallowtails along the BCS-BEC crossover. 
These results are obtained for s = 0.1 and Ep/En = 2.5. 
The quasimomentum Pcdgo = %b/2 fixes the edge of the first 
Brillouin zone. The dotted line in (b) is the half-width in 
a EEC obtained by solving the GP equation; it vanishes at 
l/kpas — 10.6. This figure is taken from Ref. [74]. 



Qualitatively, one can argue in a similar way on the ex- 
istence of the swallowtail energy band structure in Fermi 
superfluids in optical lattices; the competition between 
the external periodic potential energy and the nonlinear 
interaction energy {g\ Ui{r)v* {r)\'^) determines the en- 
ergy band structure. However, the interaction energy in 
the Fermi gas is more involved and the unified view along 
the crossover from the BCS to the BEC states are non- 
trivial and interesting problems. To answer for the ques- 
tions of 1) whether or not swallowtails exist in Fermi 
superfluids and 2) whether there exist unique features 
which make them different from those in bosons, we solve 
the BdG equations [see Eq. (9)] of a two-component un- 
polarized dilute Fermi gas subject to the one-dimensional 
(ID) optical lattice [74]. 
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FIG. 8: (Color online) (a) Lowest three Bloch bands of the 
quasiparticle energy spectrum at fcx = for P = Pedge and 
l/kpas — —0.62. Thin black dashed lines labeled by I's show 
the approximate energy bands obtained from Eq. (36) using 
fj. ~ 2.66£'r and |A| ~ |A(0)i ~ 0.54SR. (b) Bloch band of 
the quasiparticle energy spectrum around the chemical poten- 
tial for P/Pcdge = (black dotted), 0.5 (green dashed), and 
1 (red solid), respectively, at fcx = and l/kpUs = —0.62 in 
the case of s = 0.1 and Ef/Eb. = 2.5. The each horizontal 
line denotes the value of the chemical potential for the corre- 
sponding value of P. This figure is adapted from Ref. [74]. 



A. Swallowtail energy spectrum 

The energy per particle in the lowest Bloch band as 
a function of the quasimomentum P for various values 
of l/kpUs is computed [75]. The results in Fig. 7(a) 
show that the swallowtails appear above a critical value 
of l/kptts where the interaction energy is strong enough 
to dominate the lattice potential. In Fig. 7(b), the half- 
width of the swallowtails from the BCS to the BEC 
side is shown. It reaches a maximum near unitarity 
(l/fcpas = 0). In the far BCS and BEC hmits, the width 
vanishes because the system is very weakly interacting 
and the band structure tends to be sinusoidal. When 
approaching unitarity from either side, the interaction 
energy increases and can dominate over the periodic po- 
tential, which means that the system behaves more like 
a translationally invariant superfluid and the band struc- 
ture follows a quadratic dispersion terminating at a max- 



imum P larger than Pedge- 

The emergence of swallowtails in the BCS side and for 
Ep I Jc; 1 is associated with peculiar structures of the 
quasiparticle energy spectrum around the chemical po- 
tential. In the presence of a superflow moving in the z 
direction with wavevector Q (= P/h), the quasiparticle 
energies are given by the eigenvalues in Eq. (9). Since 
the potential is shallow {s <^ 1), some qualitative results 
can be obtained even ignoring Vext{z) except for its pe- 
riodicity. With this assumption we obtain 



{k, + 2qBl)Q 



n 2 



2m 



|A|2 



(36) 

with I being integers for the band index. If Q = 
0; the I = band has the energy spectrum 
\/[{k\ + k1)/2m — /i]2 + |A|2 which has a local maxi- 
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FIG. 9: Incompressibility at P = -Pedge around the crit- 
ical value of l/kpas ~ —0.62 where the swallowtail starts to 
appear. The quantity is the incompressibility of the ho- 
mogeneous free Fermi gas of the same average density. In both 
panels, we have used the values s = 0.1 and Ep/Er = 2.5. 
This figure is adapted from Ref. [74]. 



mum at kz — ~ 0. When Q = Podgo/^, the spectrum 
is tilted and the local maximum moves to fcz ~ ge/Z pro- 
vided |A| ^ E-p (and Ep/E^^ ^1)- In the absence of the 
swallowtail, the full BdG calculation indeed gives a local 
maximum at kz = ge/S and the quasiparticle spectrum 
is symmetric about this point, which reflects that the 
current is zero. As E-p /En increases, the band becomes 
flatter as a function of kz and narrower in energy. 

In Fig. 8(a), we show the quasiparticle energy spec- 
trum at A:^ = for P = Podge- When the swallowtail is 
on the edge of appearing, the top of the narrow band just 
touches the chemical potential fi [see the dotted ellipse in 
Fig. 8(a)] . Suppose 1 / kpas is slightly larger than the crit- 
ical value, so that the top of the band is slightly above fi. 
In this situation a small change of the quasimomentum P 
causes a change of ^. In fact, when P is increased from 
P = -Podge to larger values, the band is tilted and the 
top of the band moves upwards; the chemical potential /i 
should also increase to compensate for the loss of states 
available as shown in Fig. 8(b). This implies dfi/dP > 0. 
On the other hand, since the system is periodic, the ex- 
istence of a branch of stationary states with d^/dP > 
at P = Pcdgc implies the existence of another symmetric 
branch with dfji/dP < at the same point, thus suggest- 
ing the occurrence of a swallowtail structure. 



B. Incompressibility 

A direct consequence of the existence of a narrow 
band in the quasiparticle spectrum near the chemical 
potential is a strong reduction of the incompressibility 
= nd^(n)/dn close to the critical value of l/kpUs in 
the region where swallowtails start to appear in the BCS 
side (see Fig. 9). The dip of occurs in the situation 
where the top of the narrow band is just above fi for 
P — Pedge {l/kptts is slightly above the critical value). 




-0.4 -0.2 0.2 0.4 -0.4 -0.2 0.2 0.4 



z/d z/d 

FIG. 10: Profiles of the pairing field iA(2:)[ and the density 
n{z) at l/kpUs — —0.8 (blue dotted), —0.6 (green dashed), 
and —0.4 (red solid) for P = Pedge in the case of s = 0.1 and 
E-p I En — 2.5. The swallowtail starts to appear at a critical 
value of l/kpas ~ —0.62. This figure is taken from Ref. [74]. 

An increase of the density n has little effect on ^ in this 
case, because the density of states is large in this range of 
energy and the new particles can easily adjust themselves 
near the top of the band by a small increase of fj,. This 
implies that dfj,(n)/dn is small and the incompressibility 
has a pronounced dip [76] . It is worth noting that in the 
BEC side the appearance of the swallowtail is not associ- 
ated with any significant change of incompressibility. In 
fact, for the Bose gas of the average density nto of bosons, 
the exact solution of the GP equation gives k^^ = uiiqUq 
near the critical conditions for the occurrence of swal- 
lowtails, being a smooth and monotonic function of the 
interaction strength. 

C. Profiles of density and pairing fields 

Both the pairing field and the density exhibit inter- 
esting features in the range of parameters where the 
swallowtails appear. This is particularly evident at the 
Brillouin zone boundary, P = Pedge- The full profiles 
of |A(z)| and n{z) along the lattice vector (z-direction) 
are shown in Fig. 10. In general, n{z) and |A(z)| take 
maximum (minimum) values where the external poten- 
tial takes its minimum (maximum) values. By increasing 
the interaction parameter l/Zcpas, we find that the order 
parameter |A| at the maximum {z = ±d/2) of the lattice 
potential exhibits a transition from zero to nonzero values 
at the critical value of l/kpUs at which the swallowtail 
appears. Note that here we plot the absolute value of A; 
the order parameter A behaves smoothly and changes 
sign across zero. 

In Fig. 11, we show the magnitude of the pairing field 
|A(z)| and the density n{z) calculated at the minimum 
(z — 0) and at the maximum (z = ±c?/2) of the lat- 
tice potential. The figure shows that |A(c?/2)| remains 
zero in the BCS regime until the swallowtail appears at 
l/kpQs ~ —0.62. Then it increases abruptly to values 
comparable to |A(0)|, which means that the pairing field 
becomes almost uniform at P = Pedge in the presence of 
swallowtails. As regards the density, we find that the am- 
plitude of the density variation, n(0) — n{d/2), exhibits a 
pronounced maximum near the critical value of l/fcpa^. 
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FIG. 11: (Color online) Profiles of (a) the pairing field |A(z)| 
and (b) the density n{z) along the change of 1/fcpaa for P = 
Podge in the case of s = 0.1 and Ef/E-h — 2.5. The values 
of |A(z)| and 71(2) at the minimum {z — 0, blue □) and at 
the maximum {z = ±d/2, red x) of the lattice potential are 
shown. The vertical dotted lines show the critical value of 
1/fcpaa above which the swallowtail exists. The dotted curve 
in (a) shows jAj in the uniform system. This figure is taken 
from Ref. [74]. 



In contrast, in the BEC side, the order parameter and 
the density are smooth monotonia functions of the inter- 
action strength even in the region where the swallowtail 
appears. At P = Pedge, the solution of the GP equation 
for bosonic dimers gives the densities nb{0) = ni,o{l + 
Vo/2nboUo) and nb{d/2) = nbo(l - Vo/2nboUo) with 
Vo/2nboUo = i3TT/4){sEn/EF)il/kFas), where Uq = 
ATrfr^ab/nib [69, 77, 78]. Near the critical value of l/fepOs, 
unlike the BCS side, the nonuniformity just decreases 
all the way even after the swallowtail appears. The 
local density at z = d/2 \s zero until the swallow- 
tail appears in the BEC side while it is nonzero in the 
BCS side irrespective of the existence of the swallow- 
tail. The qualitative behavior of |A(z)| around the crit- 
ical point of l/fepOs is similar to that of nb{z) because 
Ubir) = (m2as/87r)|A(r)|2 [79] in the BEC limit. 



VI. QUANTUM PHASES OF COLD ATOMIC 
GASES IN OPTICAL LATTICES 

So far, we have discussed some superfluid features of 
the cold atomic gases in an optical lattice mainly focusing 
on our works. In this section we discuss some important 
topics regarding various quantum phases that the cold 
atomic gases on optical lattices can show or simulate. In 
this section, strong enough periodic potentials to ignore 
the next-nearest hopping and small average numbers of 
particles per lattice site will be assumed so that the sys- 
tem is described by the Hubbard lattice model. 



A. Quantum phase transition from superfluid to 
Mott insulator 



At a temperature of absolute zero, cold atomic gases 
in optical lattices can undergo a quantum phase tran- 
sition from superfluid to Mott insulator phases as the 
interaction strength between atoms is tuned from weakly 



to strongly interacting regime [80]. The quantum phase 
transition of an interacting boson gas in a periodic lattice 
potential can be captured by the following Bose-Hubbard 
Hamiltonian; 



H 



lb, 



h.c. 



1): 



where hi, hi are annihilation and creation operators of a 
bosonic atom on the i-th. lattice site and hi ~ blhi is the 
atomic number operator on the i-th site. The first term 
is the kinetic energy term whose strength is character- 
ized by the hopping matrix element J between adjacent 
sites and U{> 0) in the second term is the strength 
of a short-range repulsive interactions between bosonic 
atoms. 

In the limit where the kinetic energy dominates, the 
many-body ground state of N atoms and M lattice sites 
is a superposition of delocalized Bloch states with lowest 
energy (quasimomentum k = 0): 



M 



N 



I* 



SF)c//j=o(x(&Lor|0)(x^&n |0), (37) 



where |0) is the empty lattice. This state has perfect 
phase correlation between atoms on different sites while 
the number of atoms on each site is not fixed. This su- 
perfluid phase has gapless phonon excitations. 

In the limit where the interactions dominates (so- 
called, "atomic limit"), the fluctuations in the local atom 
number become energetically unfavorable and the ground 
state is made up of localized atomic wavefunctions with 
a fixed number of atoms per lattice site. The many-body 
ground state with a commensurate filling of n atoms per 
lattice site is given by 



M 



i*Mi),//c/=o«n(^^rio)- 



(38) 



i=l 



There is no phase correlation between different sites be- 
cause the energy is independent of the phases of wave- 
functions on each site. This Mott insulator state, unlike 
the superfiuid state, cannot be described by a macro- 
scopic wavefunction. The lowest excited state can be 
obtained by moving one atom from one site to another, 
which gives the energy gap of A = (J7 /2)[{n + 1)^ -I- (n — 
1)2 - 2^2] = U. 

As the ratio of the interaction term to the tunneling 
term increases ({// J can be controlled by changing the 
depth Vq of the optical lattice even without using the Fes- 
hbach resonances), the system will undergo a quantum 
phase transition from the superfluid state to the Mott 
insulator state accompanying the opening of the energy 
gap in the excitation spectrum. Greiner et al. realized ex- 
perimentally the quantum phase transition from a Bose- 
Einstein condensate of ^^Rb atoms with weak repulsive 
interactions to a Mott insulator in a three-dimensional 
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optical lattice potential [16]. Notably, they could induce 
reversible changes between these two ground states by 
changing the strength Vq of the optical lattice. The su- 
perfluid to Mott insulator transitions were also achieved 
in one- and two-dimensional cold atomic Bose gases [81]. 

Mott insulator phases of atomic Fermi gases with repul- 
sive interactions on the three-dimensional optical lattice 
have been realized and the entrance into the Mott insu- 
lating states was observed by verifying vanishing com- 
pressibility as well as by measuring the suppression of 
doubly occupied sites [82, 83]. 

B. Quantum phase transition from paramagnet to 
ant ifer r omagnet 

Sachdev et at showed that the onc-dimcnsional Mott 
insulator of spinlcss bosons in a tilted optical lattice can 
be mapped onto a quantum Ising chain [84]. The Bose- 
Hubbard Hamiltonian for a tilted optical lattice takes the 
form: 

H = -J^(S,6|+i+Sj6,+i)+^^fi,(n,-l)-£;^in, , 

i i i 

where E is the lattice potential gradient per lattice spac- 
ing. For a tilt near E = U, the energy cost moving 
an atom to its neighbor (from site i to site i -I- 1) is 
zero. If we start with a Mott insulator with a single 
atom per site, an atom can resonantly tunnel into the 
neighboring site to produce a dipole state (at the link) 
consisting of a quasihole-quasiparticle pair on nearest 
neighbor sites. Only one dipole can be created per link 
and neighboring links cannot support dipoles together. 
This nearest-neighbor constraint is the source of the ef- 
fective dipole-dipole interaction which results in a den- 
sity wave ordering. Defining a dipole creation opera- 
tor dj = bSl^j^/ \/2, the Bosc-Hubbard Hamiltonian is 
mapped onto the dipole Hamiltonian: 

H = - V2 J ^(d| + d,) + {U-E)Y, d\d^ 

i i 

with the constraint d\di < 1 and dldidl_^^di^i ~ 0. If 
the dipole present /absent link is identified with a pseu- 
dospin up/down, = — 1/2, the pseudospin-1/2 
Hamiltonian takes the form of a quantum Ising chain: 

H ^ JsJ2 StSl^i - 2V2 J + -D)Y. 

i i i 

i 

where D = E-U Siiid J2;Js{Sf + l/2)(5f+i + 1/2) 
with J5 00 is added to the Hamiltonian for imple- 
menting the constraint dldidl_^_^di+i = 0. The dimen- 
sionless transverse and longitudinal fields are defined as 
= 23/2 J/ Js h,^ = l-D/Js. 



Thereby a quantum phase transition from a param- 
agnetic phase {D < 0) to an antiferromagnetic phase 
{D > 0) can be studied by changing the lattice potential 
gradient E between adjacent sites. Using *^Rb atoms 
quantum simulation of antiferromagnetic spin chains in 
an optical lattice was made and the phase transition 
to the antiferromagnetic phase from the paramagnetic 
phase was observed [85]. 



VII. SUMMARY AND OUTLOOK 

In this article, we have focused on some superfluid 
properties of cold atomic gases in an optical lattice pe- 
riodic along one spatial direction; basic macroscopic and 
static properties, stability of the superflow, and peculiar 
energy band structure. As a complement, some phases 
other than the superfiuid phase have been discussed in 
the last section. Considering the rapid growth and inter- 
disciplinary nature of the researches on cold atomic gases 
in optical lattices, it is practically impossible to cover all 
perspectives. We conclude by mentioning a few of them. 

High controllability of the lattice potential opens many 
possibilities: Optical disordered lattice can be con- 
structed to study the problem of Anderson localization 
of matter waves and the resulting phases [86] . The time 
modulation of the optical lattice was shown to tunc the 
magnitude and the sign of the tunnel coupling in Hub- 
bard model, which enables us to study various phases 
[87]. Cold atomic gases in optical lattices might uncover 
many exotic phases which are still in debate or even lack 
solid state analogs. 

Due to their long characteristic time scales and large 
characteristic length scales, cold atomic gases are good 
playground for the experimental observation and control 
of their dynamics. Particularly, a sudden quench can be 
realized experimentally. The non-equilibrium dynamics 
after sudden quenches can be studied with high precision 
in experiments to discriminate between candidate the- 
ories. The long-standing problems such as thermaliza- 
tion, its connection with non-intcgrability and/or quan- 
tum chaos are actively sought-after topics in this direc- 
tion [88]. 

The charge neutrality of atoms seems to limit the use 
of this system as a quantum simulator. However, inter- 
nal states of an atom, together with atom-laser interac- 
tion, can be exploited for the atom to gain geometric 
Berry's phases, which amounts to generating artificial 
gauge fields interacting with these charge-neutral par- 
ticles [89] . In this way, interactions with electromagnetic 
fields, spin-orbit coupling, and even non-abelian gauge 
fields can be emulated to open the door to study the 
physics of quantum Hall effects, topological supercon- 
ductors/insulators, and high energy physics using cold 
atoms in the future. 

Cold atomic gases provide a promising platform for 
controlling dissipation in addition to engineering Hamil- 
tonian, e.g., by controlling the coupling with a subsystem 
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acting as a reservoir and by using external fields which 
induce losses of trapped atoms, etc. This possibility will 
allow us to use the cold atomic gases as quantum simu- 
lators of open systems. In addition, the controlled dis- 
sipation will offer an opportunity to study quantum dy- 
namics driven by dissipation and its steady states, non- 
equilibrium phase transitions among the steady states 
determined by the competition between the coherent 
and dissipative dynamics. Furthermore, the control- 
ling dissipation will pave a way to design Liovillian and 
dissipation-driven state preparation (e.g., Ref. [90]). 

Last but not least, cold atomic gases in optical lattices 
are also useful for the precision measurements; "Optical 
lattice clock" consists of millions of atomic clocks working 
in parallel trapped in an optical lattice. The large num- 
ber of the simultaneously-interrogated atoms greatly im- 
proves the stability of the clock and the state-of-the-art 
optical lattice clocks outperform the primary frequency 
standard of Cs clocks (Ref. [17] and references therein). 
Blocli oscillations of cold atomic gases in optical lattices 
offer a promising way of measuring forces at a spatial 
resolution of few micrometers (e.g., Ref. [91]). The pre- 



cision measurement devices/techniques will enable high 
precision tests of time and space variation of fundamental 
constants, the weak equivalence principle, and Newton's 
gravitational law in short distances. 
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